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The transport properties of a suspended carbon nanotube probed by means of a STM tip are 
investigated. A microscopic theory of the coupUng between electrons and mechanical vibrations is 
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vibron is located. This fact has profound consequences on the transport properties, which allow to 
extract information on the location and size of the vibrating portions of the nanotube. 
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I. INTRODUCTION 

Carbon nanotubes (CNTsj^ are extremely versatile 
systems with metallic or semiconducting behavior de- 
pending on their wrapping orientation.^^ Deposing them 
on an insulating substrate and tunnel-coupling it to bi- 
ased electrodes it is possible to create a single-electron 
transistor, in which the nanotube behaves as a quan- 
tum dot.^ ^ Alternately one can embed a quantum dot 
into the nanotube via geometrical defects or external 
gates, thus building a nanotube dot tunnel-coupled to 
contacts.^ 

Recent improvements in manipulation techniques have 
allowed to suspend nanotubes between two contacts. In 
this case nanotubes behave as mechanical resonators,^ 
with possible applications ranging from ultra-sensitive 
mass sensing to displacement sensors.^ Among the differ- 
ent mechanical vibration the radial breathing mode is 
the highest in energy^^ followed by the twist and the 
stretching ones. The latter have received a lot of ex- 
perimental attention Q3HIll ^Iso in view of the peculiar 
features induced on transport, such as negative differen- 
tial conductance^^ or Franck-Condon blockade. Bend- 
ing modes have usually energies lower than the experi- 
mental temperature requiring external AC drivings .'^^'^ 
Transport experiments have been employed to analyze 
the structure of nanotubes exploiting a scanning tun- 
neling microscope (STM) tip.^^ Effects such as spin- 
charge separation were observed studying the differential 
conductance as a function of the tip position. Super- 
conducting probes have been used to extract the non- 
equilibrium electron energy distribution function.!^ Also, 
the effects of chemical or magnetic impurities adsorbed 
along the nanotube were considered. 123123 
Scanning tunnel microscopy experiments have also been 
performed on suspended nanotubes. In particular, it has 
been shown how electrons injected from a tunnel micro- 
scope tip can excite, detect and control a specific vibra- 
tional mode.f^^lfi^ 

Owing to their small waist (of the order of some nm) 
nanotubes behave as a one-dimensional interacting elec- 
tronic system.^ ^ Typically, correlated quantum systems 



are studied by means of numerical techniques.'^^''^ How- 
ever, due to their inherently one-dimensional nature, car- 
bon n anotu bes are described in terms of a Luttinger 
model.'^^'^In this context, transport from a tunn eling tip 
to a static nanotube has been recently considered l^^l^^^ ^^ 
The coupling between the electrons and vibrational 
modes has been extensively studied in literature .f^^Ell 
In most cases t he si mple Anderson-Holstein model has 
been employed, in which the vibron couples only 
to the total charge neglecting the spatial modulation of 
the charge density. The Anderson-Holstein interact ion 
yields "^o^iiion- independent Franck-Condon factors^^^^ 
with visible effects in transport properties j^^tt57| 
Recently, a microscopic theory involving the coupling 
with spatial fluctuations of the nanotube electronic den- 
sity has also been developed^ in order to explain anoma- 
lous transport behaviors.!^ 



In this paper we investigate the possibility of creating 
an electrical probe for the stretching vibrational modes 
of a suspended carbon nanotube, by means of a scanning 
tunnel microscope tip. Building on the theory outlined 
in Ref. [13 we describe the coupling between vibrons and 
total charge as well as the spatial charge density modu- 
lations. This coupling gives rise to a position-dependent, 
electron- vibron coupling which strongly affects the trans- 
port properties. Position- dependent tunneling rates and 
conductance arise. This allows to obtain precise infor- 
mations about the vibrational mode of the nanotube. 
Effects are visible in metallic nanotubes and are more 
pronounced in semiconducting ones. 



The paper is structured as follows. In Sec. |II A| a Lut- 
tinger liquid model for a carbon nanotube with open 
boundary conditions is introduced. In Sees. IIBpI C the 
lattice vibrations, the electron-vibron coupling and its 



diagonalization are discussed. Section [II D| is devoted to 
the transport properties. Our results are illustrated and 
commented in Sec. lIIII Conclusions are drawn in Sec. lIVl 
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II. MODEL AND METHODS 

A. Modeling a carbon nanotube quantum dot 

The electronic properties of the CNT are characterized 
by the wrapping vector Wn,m = na-^i +ma_i, where a±i 
represent the basis vectors of the graphene latticed Due 
to the wrapped nature of the system, the energy spec- 
trum is composed of subbands corresponding to trans- 
verse excitations along the waist of the tubule. In a typ- 
ical experiment only the lowest-lying subband is occu- 
piedp^ giving a one dimensional character to the CNT. 
In this regime both n-doped semiconducting CNTs (away 
from the band gap) and metallic CNTs can be described 
in the low energy sector as Luttinger liquids with four 
branches, ^^ ^^ labeled by a = ±1, stemming from the 
two Dirac valleys of the graphene, and by s = ±1, denot- 
ing the z component of the electron spin (units h/2). 
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FIG. 1: (Color online) Schematic setup of a CNT, suspended 
between the two substrates SI and S2 at positions xi and X2. 
A quantum dot with ends at x = and x = L is embedded 
in the CNT. 



The system under investigation is schematically depicted 
in Fig. [T] it consists of a CNT, suspended between two 
substrates at Xi and X2 and free to vibrate. Embedded 
in the CNT there is a quantum dot of length L with ends 
at x = and x = L. We assume xi < and X2 > L in 
order to mimic a quantum dot inside the CNT, due to 
geometrical defect or external gates, or to treat the CNT 
itself as the quantum dot. In the following, we focus on 
the description of the quantum dot with open boundary 
conditions. 

The bosonized hamiltonian is H = Hj with 



J2^,{q)b]{q)b,{q), {h = 1) (1) 



where j G p_, cr+, cr_} are the four linear combina- 
tions of states in the a, s branches that diagonalize the 
Coulomb interaction. Here, Nj represent the zero modes 
counting the excess electrons in the j sector 

a,s a,s 

The bosonic operators bj (q) trigger collective excitations 
of the electron system with momentum q = irn/L with 




FIG. 2: (Color online) Fermi velocity vf of a semiconducting 
CNT of length L = 400 with iVo electrons and different CNT 
waist lengths |wn,m|: blue (solid) 1 nm - such as for (3,2), 
purple (dashed) 2 nm - such as for (5,4), yellow (dotted) 3 
nm - such as for (8,6), green (dash-dotted) 4 nm - such as for 
(11,7). 



n G N*. They are connected to the a^s modes by a 
Bogoljubov transformation.^^ Note that the mode j = 
p+ represents the total charge of the system. 
The collective modes propagates with velocities vj = 
Vf I gj with gp^ = g and gj = 1 Vj 7^ p+ . Here g parame- 
terizes the strength of electron scattering, with ^ < 1 for 
repulsive interactions. Note that only the velocity of the 
total charge mode is renormalized. The corresponding 
energies are ujj{q) = Vjq and Ej = nVj/AgjL. 
The value of the Fermi velocity v^ depends on the prop- 
erties of the CNT. For a metallic CNT (|n - m|/3 gN ) 
the dispersion relation is linear with v^ = vq = 8-10^ m/s 
around the degenerate point k. In the following, we will 
assume an effective n-doping with a shift of Ey towards 



higher energy values and new Fermi points /Cp 
with = E-p/v-p. 
In a semiconducting CNT the conduction and valence 
bands are separated at momentum by the direct gapP 



Attvq 

3|Wn,ri 



(2) 



with |w^^^| = aV^^ + + the waist length of the 
CNT. We will consider n-doped semiconducting nan- 
otubes with Ep > A/2, having chosen the energy ref- 
erence to lie in the middle of the band gap. Doping 
gives rise to two new Fermi points k^^ =k'^qY where 
gp 2m*V£^F - A with m* A/2vq the effective 
mass.l^ Around the new Fermi points, the Fermi veloc- 
ity is 



Vy 



3|Wn,m|A^0 

8L 



(3) 



with Aq the total number of excess electrons in the quan- 
tum dot. The dependence of v^ as a function of Aq is 
shown in Fig. [2] One observes a lower velocity with re- 
spect to the metallic case. 
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The electron field operator ^s(r) = ^^(x,?/) has to sat- 
isfy open-boundaries conditions ^^(O,?/) = ^^(L,?/) = 0. 
It can be written in the bosonized form 

= E E frA^y^" ^+hra Arx) (4) 

r=±l a^±l 

in terms of the right-movers field operators 



The coupling between electrons and vibrations can be 
microscopically derived starting fromthe tight-binding 
theory of a distorted CNT lattice.l^^ For the typical 
experimental situations one has C > 100 nm then a con- 
tinuum elastic model is appropriate with^ 



(9) 



where 



\j2i\a 



(Np^ +aNp_ +siV^^ ^asN^_ ) 



>t [0p+ {x)-\-a4)p_ (x)+s0CT^ {x)-\-as4)(7_ (a^)] 



(5) 



y^^^ 



qL 



(qx) \bj{q) + b\{q) 



+ zy^sin( 



(qx) [6,.(g)-S](q)]} 



(6) 



Here fja^s are Major ana fermions, [Oa,s^^a,s] = h and 
a is the length cutoff. The functions fr,a{^) in Eq. Q 
consist of a superposition of wavefunctions for orbitals, 
peaked around the positions of atoms in the CNT and 
oscillating with a typ ical wave vector Kq oc where 
a ^ 2.5 • 10~^^ m.'^^'^ Their specific form depends on the 
type of nanotube under consideration and will be not 
discussed here. 



Here^ c 30 eV and p(r) = ^|(r)^5(r) is the elec- 
tronic density operator, with ^s(r) given by Eq. Q. 
It consists of two components: a long wavelength part 
Plw(i*) and an oscillatory contribution psw(i*) fiuctu- 
ating on a length scale Kq^ oc a. This latter compo- 
nent, does not make sizeable contributions since qo <C 
Kq. Under the realistic assumption of strongly local- 
ized atomic orbital^^ with negligible overlapping, the 
electron- vibron coupling becomes 



Ha- 



r 



dx PLwix)d^Upix) , 



(10) 



with x< = max{0,x^o}7 ^> = min{x^i,L} and 



Plw 



1 

27r L' 



dx4'p+ (^) -\- X ^ —X 



(11) 



written here directly in its bosonized form. 



B. Lattice vibrations and electron-vibron coupling 



C. Diagonalizing the electron-vibron coupling 



We consider the case of a vibrating portion of the CNT 
(the vibron), of size £, located between x^o and Xyi. As 
confirmed in a recent experiment,^ the vibrating part 
can be different from the CNT dot. Thus, we will formu- 
late the theory in this most general case. We focus the 
description on the stretching mode. The p-th mode has 
energy cjq = wv^j L with ^ 2.4 • 10^ m/s the veloc- 
ity of the stretching modes which is approximately non- 
dispersive. In most experiments the fundamental mode 
with p = 1 is observed. 'I^'^ For these reasons, we will 
concentrate to the case of small p < 3. The p-th mode is 
described as a harmonic oscillator 



- 2 

2M 



MujI ^2 

— - — 



•> 



(7) 



where M = 27r|w^^^|>C/)o is the vibron mass with po ^ 
6.7 • 10"'' Kg/m^ the graphene density and Xq is the 
amplitude operator of the strain field 



'Up(r) = V2Xq sin 



7r(x - Xyo) 



p- 



(8) 



The latter represents a standing wave with momentum 
±qo with Qo = ptt/C. 



The relevant terms of the electron-vibron coupling are 



h = H, 



(0) 



r(pl) _ 



(Pl) 
P+ 



iJv + i^d-v with H, 



(0) 

p+ 



and i^^^'^ = Introducing = 

ibp_^{7r/jj/L) and ^J2uJp^ {T:ii/L)Xp = -\- Bj^ we have 



h = 



-E 



p2 

2M 



o 




MCoXoNp 



M>1 



(12) 



with [X^,P.] 



iS, 



jiuoi/g^ uoi = tivy/L. In 



terms of these new variables the density operator is 



271VY /7r/ix\ 

y ^ ji cos I — I 



L3 



/i>l 



(13) 



The last term in Eq. ( 12 ) describes a central harmonic os- 



cillator (vibron) linearly coupled to a infinity of harmonic 
oscillators (plasmon modes of the dot). Note that, for 
reasonable experimental parameters and considering the 
lowest stretching modes, one always has uo^ > uq both 
for metallic and semiconducting CNTs. Additionally, the 
vibron is also coupled to the total average charge Np^ on 
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the quantum dot in analogy to the Anderson-Holstein 
model. The coupling coefficients are 

Co = V2XmUjl^'^Jo ; C^>i = 2XmUjQ^'^y^J^ (14) 
where {hz > 0) 



with 



J. 
and 



dx cos 



Am — 



tvTTX 



L J 



cos 



^S\/P07r|wn,m|^s 



{X - Xyo) , (15) 



(16) 



Taking as a reference^ ^ a waist length of about 2 nm, 
one has Am ^ 2. Note that a recent experiments^ reports 
a larger c which leads to a larger Am- One finds 

J, = j(°) + 4-)0(-O+-x„o) + 4+)^(xi-L-O+), (17) 

with 9{x) the Heavyside step function and 



4°^ 



7r(p^ — K^5'^) 

pS sin (/rp^o /^) ~ sin (tt/t^^q) 

7r(p^ — tz'^S'^) 
{-1)Pk6'^ sin [7rA>:(^o + ^^)] 

7r(p^ — tv'^S'^) 
pSi-l^sin [7TpS{l-^o)/S] 



where 5 = C/L and = Xy^/L. 

Let us now comment the general features of the cou- 
pling depending on the relative size and position of dot 
and vibron. When the vibron is much larger than the 
dot {5^1) and the latter is embedded into it one has 
Jk, ^ namely in the large vibron limit the coupling to 
the charge density fluctuations vanishes and only the con- 
ventional Anderson-Holstein coupling survives. In this 
limit, the electron- vibron coupling reduces to the stan- 
dard fornPlcJoAm^o ^-^o^'p+ with = ^/Muj^, and the 
coupling constant is Am- 

The most interesting case occurs when the vibron is 
smaller than the dot and embedded into it. In this 
regime, one finds Jq = 0, while J^>i ^ 0. This fact 
signals a radical departure from the Anderson-Holstein 
model: the coupling between electrons and vibrons oc- 
curs only via the spatial fluctuations of the electron den- 
sity. 

We now turn to the diagonalization of the electron- vibron 
coupling. The terms linear in in Eq. (12) can be 
exactly diagonalized,'^ leading to 




(18) 



z/>0 



(19) 



In order to diagonalize the term (x TVp^ a Lang-Firsov 
canonical transformation is used 

^ = g-i^P+E.>o^^^^^ (20) 

with r]y = \fMCQkoj,/(Vtl\fM). This leads to a shift 

X,^X,- rj^N,^ , (21) 
which finally casts the hamiltonian into the diagonal form 



/i>0 



2 



with ^E,^=ClY..^,{klJ^ll). 

The energies Vt^ of the new eigenmodes are the roots of 
the secular equation 



= ^0 + E 



(23) 



V>1 



with 



h 



-^^J±^ko^ (withM>l), (24) 



i + E 



r^2 



-1/2 



(25) 



As can be clearly seen, the modes with /i > 1 have 

1.02r 




FIG. 3: Energy of the lowest eigenmode. (a) Plot of Q.q/ujq 
as a function of a; (b) Plot of Q.i/uJi as a function of a. In 
all figures, ^ = 1, Xvo — 0^ p — 1^ g — 1 and Am — 2. 



energies Vt^ > jiuoi and represent blue-shifted dressed 
plasmons. The lowest-lying solution, on the other 
hand, has VLq < uq and represents a dressed vibronic 
mode red-shifted by the electron-vibron interaction. By 
inspecting Eq. (23) one always obtains a real solution 
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for Qjj^ with /i > 1. On the other hand, the existence 
of a real solution for fto requires oOqUjI > X]/2>i C'^/m^- 
When this condition is not fulfilled, the Wentzel-Bardeen 
instability occurs. In our calculations we have always 
checked that for realistic parameters the system does 
not exhibit this instability. The energy of the dressed 
vibronic mode is very sensitive to the ratio a = v-p/vs 
between the Fermi and the sound velocity. While for 
a = 32, corresponding to the case of a metallic CNT, 
one has fto ^ ujq, for lower values of a, typical of a 
semiconducting CNT, a suppression of occurs, see 
Fig. [3]^ a). Note that the dressed plasmons are almost 
insensitive to the ratio a (cf. Fig. [sj^b)). 



The transformations in Eq. jTol ) and ( [2Q| ) affects the elec- 
tronic field operator of Eq. (|5|. Up to an irrelevant phase 
constant, the field (pp^ix) is 



/x>0 



where 



= V2cji 2^ Kp sm \ -^] 

V cji ^-^ V \ 



2 X ^ 

D. Modeling transport 
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FIG. 4: Schematic setup of a STM transport experiment per- 
formed on a suspended CNT. The quantum dot is the nan- 
otube itself, biased with respect to the tip and two lateral 
contacts. Electrons can flow through the system via the STM 
tip at position (x) and at the contacts (xi,2). A back gate 
allows to tune the effective charge on the dot. 



Our task is to model an STM tunneling tip on a 
CNT.EJtl^ The setup is sketched in Fig. [i] and is com- 
posed of a STM tip and two lateral contacts tunnel- 
coupled to a suspended CNT. Tip and contacts are biased 
in such a way that, for V > 0, electrons flow from the tip 



to the contacts through the CNT. Additionally, a back 
gate is capacitively coupled to tune the effective charge 
on the CNT-dot. From now on, we will focus on the 
most interesting regime, namely that of a quantum dot 
along all the CNT with a vibron imbedded into it. The 
hamiltonian for the vibrating CNT is then 



^CNT 2^P^ 



E 

/i>0 



p2 



2 ^ 



where A^g represents the charge induced by the back gate 
voltage V^. The dot is laterally coupled to the two Fermi 
contacts, via the tunneling hamiltonian'^ 



j — 1,2 a,s,q 



,)c,-,(g)+h.c., (29) 



where to is the tunneling amplitude, Cj^s{Q) are the op- 
erators for an electron with momentum q and spin s in 
the non-interacting lead j and xi = 0^ X2 = L are the 
position of the tunneling contacts. 

The STM tip is modeled as a semi-infinite Fermi con- 
tact, placed above the CNT at a position x along it. The 
tunnel coupling is expressed in terms of the Fermi field 
operator for the forward modes of the tipl ^^ l ^^ l '4^s,f{z) {z 
is the coordinate along the tip with z = at the vertex) 



T,T- 



/ d^ ^o(^-^)Xl^+i,«,.(^) 



V^.,F(0+)+h.c.. 



The function ro{x) = ro(f{x) describes the geometry of 
the tip, with < (p{x) < 1 peaked around x = 0, where 
(p{0) = 1. In the following we will assume a typical tip, 
with an effective width of a few atomic cells of the CNT. 
A voltage —V/2 is applied to the STM tip, while the lat- 
eral contacts are kept at the same voltage V/2. Voltage 
drops are assumed to occur symmetrically on the CNT. 
More general potential distributions do not affect the re- 
sults at a qualitative level. 

We will consider the sequential tunneling regime, treat- 
ing the tunnel couplings to the lowest perturbative order. 
Since we are interested into the low-energy transport 
regime {e\V\kBT ^ Qq) we disregard the dynamics of 
the modes /i > 1. Furthermore, we will consider the 
relevant situation of a damped vibronic mode, with a 
thermal equilibrium distribution at temperature T. 

The eigenstates of the suspended CNT can be expressed 
by \{Na^s}) specifying the distribution of excess elec- 
trons in the channel a with spin s. We will consider 
the resonance between the state of a closed shell with 
A^o = 4/^, {hz integer) electrons and zero excess charges, 
denoted as |0) = |0, 0,0,0), and Nq -\- 1 electrons in the 
state Note that no qualitative difference in our re- 
sults would occur, for a different value of A'o- There 
are four states with A'o + 1, electrons, denoted by |a, 5), 
each with one extra electron in the state a, s. They are 
all degenerate, with energy Ep_^{l — A/g)^/2 + ScJi/S. 
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We set up a master equation for the reduced density 
matrix, obtained tracing out the leads and vibron de- 
grees of freedo m - thus neglecting coherences among vi- 
brational states!^ ^^ * ^^ * Upon the assumption of a STM 
tip width of some unit cells, coherence effects between 
states \a,s) and \a^^s) {a ^ a^) are vanishing. Co- 
herence between different spin states is also absent in 
view of the absence of spin correlations in the contacts. 
Therefore, the master equation reduces to a standard rate 
equation for the occupation probability of the quantum 

dot PN,{t) = P|0)(t); PiVo+lW = Ea,sP\a,s){t)- The 

steady-state current is then written, to lowest order, in 
terms of tunneling rates 



p(C)-p(T) 



(x) 



r(c)p(T) 

J- J- /->iif 



(x) 



rin(x) + rout(^) 



with 



rin/out(^) — r 



in/out 



(x) 



.(C) 



i/out 



(30) 



(31) 



and 



(x) 



■ l«,s)^|0) 



r(c) 

out — 



the rate associated to the tip (T) and to the contacts (C). 
We quote here explicitly the expressions for the tunnel-in 
processes, the tunnel-out rates are similar. 

The lateral contact and tip rates are respectively 



|o>^l«,^> = ^oY.Biixj)f{AE + mo + eV/2),i32) 

l>0 



r 

\0)^\a,s) 



{x) = T^^^Y.Bi{x)f {AE^ mo- eV/2) {33) 



l>0 



Here, Fq = 27ri^o|^op7 Tq^^ = 27ri^okop7 is the leads 
density of state and f{E) the Fermi function. The rates 
are then a superposition of Fermi functions at energies 



p+ 



-UJl , 



(34) 



shifted by the energy IQq^ representing the contribution 
of a transport channel exciting / vibron quanta. The 
weights Bi{x) considered in the regime /cb^ <C ard^ 



Bi{x) 



X^^ix) 



where 



X\x) 



2^0 



al{x) 



^0 q2 



(35) 



(36) 



represents the l ocal electron-vibron coupling strength^^^ 
see Eqns. ( 27|28 ). This is in sharp contrast to the case 



of the Anderson-Holstein mo del ,'^^'^^1^^ appropriate for 
large vibrons. In that case, as discussed above, cou- 
pling to the density fluctuations would vanish leading 
to ao{x) = Po{x) = and to a coupling simply given by 
Amax which is clearly independent of the tunneling posi- 
tion. 



III. RESULTS 
A. Local electron-vibron coupling 

Let us analyze in details the space dependence of the 
electron-vibron coupling X{x) which determines the cur- 
rent behavior. The coupling X{x) depends both on ge- 
ometrical parameters Xyo (the vibron origin), S = C/L 
(vibron length) and on physical ones a = vf/v^^ Am de- 
fined in Eq. ( 16 ) and the electronic interaction parameter 
9- 



We remind that the parameter a is affected by the metal- 
lic or semiconducting nature of the CNT. In a metallic 
CNT one finds a = 32. On the other hand, for a semicon- 
ducting CNT lower values of a are possible (see later). 
Figure [s] shows X{x) for different vibron configurations 



0.8 



(a) 



0.6 



(b) 




FIG. 5: (Color online) Local electron-vibron coupling X{x) as 
a function of the tip position x for: (a) (5 = 1, Xvo = 0, p = 1; 
(b) as above but p = 2; (c) ^ = 0.4, Xvo = 0.15L, p = 1; 
(d) as above but p = 2. In all panels red (solid) lines denote 
a metallic CNT with a = 32, blue (dashed) lines denote a 
semiconducting CNT with a = 5. Other parameters: Am = 2 
and g = 1. The shaded plots on top of the panels depict the 



amplitude of the strain field Up 



The arrow in Panel fc) 



denotes the tip position for the conductance shown in Fig. [S] 



and CNT types. It can be seen that the electron-vibron 
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coupling strength for a metallic CNT (solid red lines) is 
smaller than that for a semiconducting CNT. Indeed, for 
a semiconducting CNT, the velocities of the electronic 
and vibronic subsystems are closer, which implies a more 
favorable interplay between them. In the rest of the pa- 
per, we will choose a = 5 to model a semiconducting 
CNT and a = 32 for the metallic one. 
The amplitude of the electron-vibron coupling is maxi- 
mal in the region where the strain field is maximum. In- 
deed, X{x) closely follows the amplitude of Up{x), which 
is sketched on top of the panels of Fig. |5] For J < 1, 
this implies a particularly sizeable X{x) only in the re- 
gion where the vibron sits, see Figs. |5jc,d). Coupling to 
higher vibronic modes produces more oscillations, as can 
be seen in Figs. [5|b,d). It also makes the electron-vibron 
coupling strength weaker. Figure [6] shows the comparison 
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FIG. 6: (Color online) Local electron vibron coupling A(x) as 
a function of x for 5 = 1, x^o = and different vibron modes: 
red (solid) p = 1, green (dashed) p = 2, blue (dotted) p = 3. 
The thin lines are the maxima of A(x). Other parameters: 
a = 5, Am = 2, ^ = 1. 



neling rates. Figure [t] shows the ratio Al/Ar as a func- 
tion of 5 for a vibron with origin at XyQ = O.IL. At small 
values of 5 the vibron is asymmetrically located near the 
left tunnel barrier. As a consequence, Al > Ar. This 
mechanism is at the origin of the systematic suppres- 
sion of conductance traces in a recent experiment. For 
increasing 5 the situation evolves towards a more sym- 
metric setup and indeed for 5 = 0.8, corresponding to a 
symmetric vibron with respect to the CNT, one recovers 
Al = Ar. 



B. Transport properties 

Figure [S] shows the density plot of the differential 
conductance Q = dl/dV in the (V^Ng) plane for the 
situation depicted in Fig. [sjc), semiconducting case. 
The large white areas at small V are the Coulomb 




between the first three vibronic modes. The intensity of 
the electron-vibron coupling strength decreases with the 
increasing order of the vibronic mode. Denoting Xp the 
maximum of X{x) for the p mode, we find Xp = Xi/^/p. 
We now briefiy comment on the value of the electron- 



FIG. 8: (Color online) Semiconducting CNT. Differential con- 
ductance Q (units e^FQ^'^/Qo) as a function of A^g and V (units 
^o/e) for a vibron originating at Xyo = 0.15L with S = 0.4, 
p = 1 and a tip at x = 0.35L. Other parameters: Am = 2, 
a = 5, ^ = 1, /cbT = 0.15 Qo and A = FoM^^ = 100. 




blockade regions where transport is interdicted and the 
CNT is occupied by Nq or A/q + 1 electrons. Within 
the transport region, delimited by the two most intense 
conductance traces, a series of equally spaced lines are 
clearly visible. They correspond to the excitation of 
the vibronic mode at energy Qq. We consider different 
tunneling amplitudes through the tip and the lateral 
contacts, introducing the asymmetry A = Fq/Fq^^ 



FIG. 7: (Color online) Ratio of the coupling strengths Al/Ar 
(see text) as a function of 6 for a vibron with origin at 
Xvo = O.IL and p = 1. The shaded plots schematically depict 
the amplitude of the strain field of the vibronic mode Up{x). 
Other parameters: a = 5, Am = 2 and g = I. 



We will concentrate the discussion on the regime k^T < 
r^o, quoting for simplicity analytical expressions for A ^ 
1 only, realistic in an STM experiment. Exploiting the 



fact that Al,r <C 1 we assume Bi{xj) ^ ^^,o 
regime {V 0) the conductance is then 



In the linear 



vibron coupling at the position of the tunneling barriers 
Al = A(0) and Ar = A(L), which govern the lateral tun- 



Ain 



2pe^T^^^Bo{x) 



cosh 



(3AE 
2 



iln(4) 



cosh 



/3AE 



(37) 
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with AE in Eq. (|34|. The logarithmic factor ln(4) stems 
from the fourfold degeneracy of the state with A^o + 1 
electrons. The amplitude of the linear conductance is 
modulated by the factor 



Bo{x) 



(38) 



Therefore, the conductance is suppressed in the region 
where the electron- vibron coupling is large. 
In the nonlinear regime (V > /cbT) one finds 



^nonli] 



2 ;>o cosh 



Biix) 



AE+ino-eV/2 
2 



(39) 



Equation ( |39[ ) represents a fan of equally-spaced conduc- 
tance peak lines located at Ep^{l — 2N^) + 3a;i/8^p^ + 
IVLq — eV/2 = 0, thus with negative slope in the (V, A/g) 
plane, see Fig. [8] They originate by the tunneling from 
the STM tip to the CNT, triggering vibronic excitations. 
We note that, because of the smallness of Al,r, the trig- 
gering process due to the tunneling on the contacts bar- 
riers is strongly suppressed (since Bi{xj) ^ di^o) with a 
corresponding absence of conductance lines with positive 
slope. Each of the above peaks is weighted by Bi{x) 
which in turn conveys informations on X{x). This is par- 
ticularly clear for the l-th {I > I) nonlinear conductance 
peak. Indeed, for X{x) < 1 (which is the case considered 
in our calculations, see Fig. [5| one finds that Bi{x) is a 
monotonically increasing function of X{x) - see Eq. ( [35| ) 
- which implies an increase of the conductance for in- 
creasing coupling strength. This fact allows to directly 
map the spatial modulations of the nonlinear differential 
conductance into modulations of X{x). 
To be more specific, let us consider the resonance case 



AE = 



TV. 



36^1 

8^^ 



(40) 



p+ 



and study Q. As is clear by inspecting Figs. [9]^a-d), the 
differential conductance Q exhibits position-dependent 
modulations in close agreement to the behavior of X{x). 
The most striking features show up indeed near the vi- 
bron position, where Q is suppressed for V ^ 0, and is 
enhanced for eV ~ 2lQo (/ = 1,2,...), see Fig. [s] 
For the case of a metallic CNT, shown in Fig.[lOj the spa- 
tial modulations of the conductance are less pronounced 
due to the decreased intensity of the electron- vibron cou- 
pling with respect to the semiconducting case. 
The close resemblance of the spatial modulations of the 
non linear Q with X{x) is supported studying the conduc- 
tance at fixed bias shown in Fig. [TT] for eV = 2Qo and 
eV = 41^0 • Clearly Q is enhanced where X{x) is large. 
This confirms that position-resolved conductance maps 
are source of valuable informations about the intensity 
of the strain field along the CNT and consequently on 
the location and size of the vibron mode. 




FIG. 9: (Color online) Semiconducting CNT: Plot of ^ (units 
e'^r^^-' /Qo) as a function of the tip position x and bias voltage 
V (units ^o/e) at resonance A^g = (1/2) + (SuJi/SEp^) and 
(a) S = 1, Xvo = and p = 1; (b) same as in (a) but for 
p = 2; (c) 5 = 0.4, Xvo = 0.15L and p = 1; (d) same as in 
(c) but for p = 2. Other parameters: Am = 2, a = 5, ^ = 1, 
A:bT = 0.15Qo and A = 100. 




FIG. 10: (Color online) Metallic CNT: Plot of G (units 
c^Fq^'^/Qo) as a function of the tip position x and bias voltage 
V (units ^o/e) at resonance A^g = (1/2) + (SuJi/SEp^) and 
(a) S = 1, Xvo = and p = 1; (b) same as in (a) but for 
p = 2; (c) ^ = 0.4, xo = 0.15L and p = 1; (d) same as in (c) 
but for p = 2. Other parameters: Am = 2, a = 32, g — 1^ 
kBT = 0.15 Qo and A = 100. 
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X 1 



FIG. 11 



(Color online) Differential conductance Q (units 
c^V'q'' as a function of the tip position x for A^g = 
(l/2) + (3a;i/8£;p+) and eV = (red solid line) or eV = 
(blue dashed line) for a metallic CNT: (a) 5 — 1, Xvo — 
and p = 1; (b) same as in (a) but for p = 2; (c) ^ = 0.4, 
Xvo — 0.15L and p — 1] (d) same as in (c) but for p — 2. 
Other parameters: Am = 2, a = 32, ^ = 1, k^T — 0.15 ^^o 
and A 100. 



C. Interaction and asymmetry effects 

We close by briefly commenting about the role of 
electron-electron interactions and that of the barriers 
asymmetry 




FIG. 12: (Color online) (a) Conductance Q (units e^r[,^V^o) 
as a function of the tip position x for decreasing values of g 
from g — 1 (noninteracting, lightest gray) to ^ = 0.2 (strongly 
interacting, black); (b) Plot of Q for different asymmetries (see 
key) and g — 1. Other parameters: ^ = 1, x^o = 0, p = 1, 
Am = 2, a = 5, /cbT 0.15 Oq and eV 2Qo. 



Figure [T2|a) shows the position-resolved nonlinear con- 
ductance at eV = 2Qo for increasing Coulomb interaction 
strength (decreasing values of g). Clearly, conductance 
is suppressed in turns, signaling the suppression of X{x). 
This fact can be explained in terms of an increase of the 
velocity of the charged mode Vp^ = vy / g which induces 
an effective parameter a higher than that of the nonin- 
teracting case. 

Concerning the role of the asymmetry between the 
tip rate and the leads rate. Fig. [l2|b) shows that the 
amplitude of the conductance modulations is decreased 
when making the tunnel barriers more symmetric with 
a collapse when reversing the asymmetry A < 1 (i.e. 
making the contacts more opaque than the STM tip). 
This can be understood by observing that for A ^ 1 
the current is dominated by the slowest barrier which 
is the STM one. This implies that the space-dependent 
tunneling rate can be efficiently probed in this regime. 



IV. CONCLUSIONS 

In this paper we have shown how transport measure- 
ments, performed with a scanning tunnel microscope tip 
on a suspended carbon nanotube, can bring information 
about its vibrational stretching dynamics. This theory 
predicts a position- dependent coupling constant, which is 
larger in the region where the vibron is located. 
The position-dependent coupling constant strongly af- 
fects the tunneling rate of electrons through the tip and 
has relevant consequences in the transport spectra. In 
particular, we showed that conductance maps in the lin- 
ear and nonlinear regime, obtained sweeping the tip along 
the nanotube, are closely connected to the local coupling 
constant and allow to localize the position and size of the 
vibron. Effects can be more pronounced in semiconduct- 
ing nanotubes due to the reduced Fermi velocity which 
matches more closely the speed of the vibrational mode. 
The role of electronic interactions and of the asymmetry 
between tip and metal contact tunnel barriers have also 
been addressed. 

This work could inspire a new class of experiments which 
aim at studying the vibrational degrees of freedom of a 
nanotube by means of electrical measurements. 
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